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^ ■ Abstract 
<3\ 

\^ ' scheme based on Ashtekar variables. A comparison with the covariant derivation 

■ and various subtleties in the calculation and choice of gauges are pointed out. Never- 

O , theless, the treatment is more systematic when correction terms of canonical quantum 

gravity are to be included. This is done throughout the paper for one example of 



Cosmological perturbation equations are derived systematically in a canonical 



^J^' characteristic modifications expected from loop quantum gravity. 

ix^ ■ 1 Introduction 

The backbone of most of current cosmology is the theory of perturbation equations for 
metric modes around an isotropic space-time p. It is used in particular for cosmological 
structure formation and for testing alternative theories beyond general relativity such as 
quantum gravity candidates. The underlying equations of typical interest are the linearized 
Einstein's equations, and so it is straightforward to include corrections if they come from a 
Lagrangian modified by quantum or other effects. This is, in fact, the situation encountered 
in most studies of so-called trans-Planckian issues for the effect of quantum gravity on 
structure formation. Modifications derived from a Hamiltonian formulation as it is used in 
canonical quantizations can, however, not be implemented in this direct way. Since several 
effective modifications to Hamiltonians have been derived in recent years in particular 
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within the framework of loop quantum gravity [3 El IH IS]; it is of interest to re-derive 
cosmological perturbation equations in a purely Hamiltonian fashion starting from the 
gravitational Hamiltonian. This is done in detail in this paper in a derivation based, as 
loop quantum gravity itself, on real Ashtekar variables fHlIZI- We will present a detailed 
derivation for scalar modes in longitudinal gauge around a spatially flat model, pointing out 
several subtleties compared to the Lagrangian derivation. Our analysis treats gravitational 
and matter terms on the same footing, showing how all of them can be obtained from 
the total Hamiltonian. This presents the classical basis for a systematic investigation of 
effective perturbation equations and cosmological implications resulting from canonical 
quantum gravity. In our calculations, only one type of corrections (from inverse powers of 
metric components in the Hamiltonian) is used, and their implications are discussed. 

2 Variables and Equations 

To set up linear metric perturbations f^, one perturbs the background metric 

ds^ = a\r]) (-dr/2 + ^^^dxMx'') , (1) 

here chosen as a fiat isotropic metric written in conformal time rj and with spatial coordi- 
nates x"". There are initially ten perturbation functions for the ten metric components, but 
some of them can be absorbed simply by redefining coordinates. The remaining functions, 
in gauge-invariant combinations, comprise scalar, vector and tensor modes. We are here 
primarily interested in scalar modes which in longitudinal gauge lead to a perturbed metric 

ds^ = a\r]) (-(1 + 2(l))dr]^ + (1 - 2^)6abdx''dx'') (2) 

which is thus diagonal. Moreover, in the absence of anisotropic stress it is consistent with 
longitudinal gauge to set cf) = ip, reducing the perturbations to a single function. We will 
also do so in our final equations, but not immediately since it is a consequence of equations 
of motion and should not be used in the process of deriving such equations. Expanding 
Einstein's equations to linear order in ip then leads to 



= -V5r° (3) 
a 2 



^ + 3_^ + 2 + = -a'bTl (4) 

a \a J 2 

+ = -^a^^ro (5) 

where k = SvrG is the gravitational constant and a dot denotes a derivative by conformal 
time 7]. The source terms on the right hand side of these equations are components of 
the energy momentum tensor provided by the matter ingredients, also perturbed linearly. 
These components follow from functional derivatives of the matter Lagrangian by metric 
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components, for a scalar field with potential V{ip) and Lagrangian 



we have, for instance/ 



^ (^ip6^ - ^'^ ip + a'^V,^{<^)6ip^ (6) 

ST^ = -\^5cp,a (7) 
ST,^ = 1 (-^ 50 + ^2 ^ + a' V,^ (y.) 5^ . (8) 

In addition, the matter Lagrangian determines equations of motion for matter fields 
such as those of a scalar ip: 

ip + 2^^ + V^i^) = (9) 
is the background Klein-Gordon equation, whereas 

6^ + 2-6ip - V^6ip + aV,^^ {ip)6ip + 2aV,^ ((^)V' - 4(^V' = (10) 
a 

describes the perturbed part of the scalar field. 
2.1 Canonical formalism 

Cosmological perturbation equations are the Einstein's equations expanded in metric per- 
turbations. Once a gauge is chosen and modes of interest are selected, the perturbed metric 
is specified and ready to be inserted in the expansion. Since the canonical formalism is 
equivalent to the Lagrangian one which yields Einstein's equations as the Euler-Lagrange 
equations of the Einstein-Hilbert action, the same perturbation equations must result. 
However, some of the derivations are more subtle since one has to fix gauges and select 
modes at the right places. Moreover, one first starts with a different set of variables and 
first order differential equations, which are combined to the usual second order equations. 
Keeping in mind that quantum gravity can lead to several modifications it is helpful to go 
through the canonical derivation in detailed steps, which is what we do in this paper. 

In a canonical formulation jS], the Hamiltonian H rather than Lagrangian L is the basic 
dynamical object, determining equations of motion of any phase space function / by means 
of Poisson brackets, / = {f,H}. The Poisson structure defines the kinematical arena, 
which is usually written in terms of a set of basic canonical variables such as position and 
momentum in mechanics. While dynamics as well as expressions for momenta follow from 
the same object in a Lagrangian formulation, they are separate in a Hamiltonian one. The 



^Note that the sign of the relation between pressure and T° depends on the signature of the metric. 
Our signs correspond to positive spatial metric components as in |j2Jl. 
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Poisson structure is thus prescribed independently of the Hamiltonian, but both of them are 
needed to determine dynamics. Basic configuration variables in a Lagrangian formulation of 
gravity are the components of the space-time metric Qab, and their momenta are determined 
as usually by derivatives vr"* = 6L/6gab- The dot refers here to a time coordinate in which 
the action is written. Since general relativity is covariant under arbitrary changes of 
space-time coordinates, the choice of time does not play a physical role. Nevertheless, 
by definition of its kinematical objects a canonical formulation does not appear manifestly 
covariant. Indeed, not all components of the space-time metric appear equal: some of them, 
the time-time component N and the time-space components A^*^ do not occur as first order 
derivatives in the action such that their momenta vanish identically, ttn = = TT^a. 

This is a consequence of general covariance and implies the existence of constraints. 
Since momenta of and A^"^ vanish, their equations of motion imply 



as constraint equations on the remaining phase space variables. In fact, because there is no 
absolute meaning to the time coordinate at all, the total Hamiltonian is a sum of constraints 
H = H[N] + DIN"] with the Hamihonian constraint H[N] = J d^xN{x)6H/6N and 
the diffeomorphism constraint D[N"-] = J d^xN"'{x)SH/SN°'. Coordinate time evolution 
through Hamiltonian equations of motion is completely specified only when N and A^" 
are known as functions on space-time. However, there are no equations of motion for A^ 
and A^" themselves; they are not dynamical since their momenta vanish. They have to be 
chosen in order to fix the gauge in which space-time properties are computed. That the 
constraints generate coordinate changes can more easily be seen for the diffeomorphism 
constraint which satisfies {f,D[N"-]} = C^af for any phase space function /, where on 
the right hand side the Lie derivative occurs. Both constraints receive contributions from 
gravitational fields (the spatial metric and their momenta) and matter fields. 

The Hamiltonian formulation is thus based on phase space coordinates given by the 
spatial metric components, matter fields and their momenta. In addition, there are the 
lapse function A^ and shift vector A^*^ which need to be chosen for a particular gauge. Their 
dynamical behavior is given by Hamiltonian equations of motion derived through Poisson 
brackets with the constraints. Since the Hamiltonian is usually, and in particular in gravity, 
a quadratic polynomial of the momenta conjugate to metric components, Poisson brackets 
between configuration variables and the Hamiltonian are linear in momenta. Thus, the 
Hamiltonian equations of motion for configuration variables relate momenta to first order 
time derivatives. Equations of motion for the momenta can then be reformulated as second 
order differential equations for configuration variables which agree with the Euler-Lagrange 
equations. Moreover, one can replace momenta in the constraint equations by first order 
derivatives of configuration variables, giving additional first order differential equations. 
The set of equations ©, O and (0) thus consists, from the Hamiltonian perspective, of 
two constraint equations, the Hamiltonian constraint Q and the diffeomorphism constraint 
(0) and one equation of motion (j^ for the single scalar mode. 











ttn = {vtat, H} = 

TT^a = {lT]\ra , H} 



6H/6N and 
-6H/5N'' 
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However, this set of equations is not the most general one for a hnearized metric. Gauge 
choices and a selection of modes have been made, the latter excluding vector and tensor 
modes and equating the lapse perturbation with the scalar mode. These put conditions 
on the variables and on the multipliers and A^". In a Hamiltonian formulation one has 
to be careful about when to make such choices in the process of deriving the equations 
of motion. Gauge choices have to be made from the start because this determines what 
the time variable and other coordinates in the resulting differential equations mean. For 
instance, the homogeneous mode of the lapse function has to be specified, which is usually 
chosen as iV = 1 for proper time or iV = a for conformal time. But this is to be done 
only for equations of motion, not for a derivation of the constraint. It is clear that setting 

= 1 or equal to the scale factor does not result in the right Hamiltonian constraint 
6H/6N which requires to be an independent variable. Similarly, one often sets the shift 
vector to zero, while the diffeomorphism constraint SH/SN"' is to be imposed fully. 

The correct procedure is as follows: To derive constraint equations, no gauge choices 
are to be made. In fact, without knowing the constraints it is impossible to know what 
gauge freedom one has. In the next step, one derives Hamiltonian equations of motion for 
the phase space variables, qat and their momenta in the case of gravity. Here, the gauge has 
to be chosen before computing Poisson brackets to give meaning to coordinate derivatives. 
One can also restrict to specific modes, but other conditions are not to be done. For 
instance, equating the lapse perturbation to the scalar mode is only justified as the result 
of equations of motion. Doing this before computing Poisson brackets would introduce 
erroneous relations between independent degrees of freedom. Thus, such a simplification 
must be made only in the final expressions for equations of motion. 

2.2 Perturbed canonical variables 

Also the set of canonical variables matters for a quantization: while classically one can 
change variables by canonical transformations, their quantum representations can appear 
very different. Loop quantum gravity crucially depends on properties of Ashtekar variables 
inilZj due to their transformation properties. First, one introduces a co-triad instead of 
the spatial metric qab, related to it by e^^el = gab- (Unlike the position of spatial indices 
a,b,..., the upper or lower positions of indices i are not relevant, and summing over i 
is understood even though it appears twice in the same position.) An oriented co-triad 
contains the same information as a metric but has more components as it is not a symmetric 
tensor. This corresponds to freedom one has in rotating the triple of triad co-vectors which 
does not change the metric. Not being of geometrical relevance, this freedom is removed in a 
canonical formalism by implementing the Gauss constraint introduced below. By inverting 
the matrix [el), one obtains the triad e^, a set of vector fields related to the inverse metric 
by e"e^ = q"'^. Just as the metric determines a compatible Christoffel connection P^^, a 
triad determines a compatible spin connection 

P^ = _,«.^e$(V5 + ie^eL9[eel]). (11) 
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Its components define the Ashtekar connection = + 'jK^^ together with those of 
extrinsic curvature 

= e^Ka, = -^eliCtqat + D^^N,)) (12) 

where the right hand side uses lapse function and shift vector A^*^ in addition to the 
spatial metric qab- Moreover, Ct denotes a Lie derivative along a timelike vector field chosen 
to describe changes in coordinate time, and Da is the covariant derivative compatible with 
the spatial metric Qab- In A^, we have the positive real Barbero-Immirzi parameter 7 [Hin] 
which we keep here for generality although it will not play a large role later on. 

The Ashtekar connection is thus a measure for curvature, and spatial metric information 
is described by the densitized triad E°- = |det e^| ef obtained from the triad. By multiplying 
the triad by the determinant of the co-triad, which is identical to the square root of 
the determinant of the spatial metric, it becomes canonically conjugate to the Ashtekar 
connection: 

{A^aix),E'^iy)} = ^K6l6}6ix,y). (13) 

This follows from the gravitational action which in a first order formulation contains time 
derivatives of the connection in the term 

1. [ d'x^E^ 
-fn J dt 

showing that the connection and the densitized triad are canonically conjugate. 

An unperturbed isotropic triad and connection for flat spatial slices [lOj can always 
be chosen of diagonal form E^ = p6^ and = c5\ with canonically conjugate c and p: 
{c, p} = |7/t. In more familiar terms, these variables are related to the scale factor a by 
IpI = a? and c = 7da/dt. With inhomogeneous perturbations, the triad and connection 
components will become space-dependent, but not in a way that is completely unrelated 
between E^ and A\ because of the spin connection. This implies that not both the triad and 
the connection can remain diagonal even when only scalar perturbations in longitudinal 
gauge are considered for which the spatial metric is diagonal. (This happens generally 
in inhomogeneous situations; see also fll.) Another way to see this is by looking at the 
Gauss constraint 

daE^ + e'^^E^jA^ = (14) 

which ensures invariance of physical results under rotations of the triad, which do not 
change the metric. Were both E^ and A\ diagonal, the second term e^^^EjA'^ would 
vanish, constraining inhomogeneity by daEf = 0. 

It is most useful to keep the triad diagonal since this simplifies the classical calculations, 
and even more so the quantum ones where currently only situations of diagonal triads are 
sufficiently accessible by explicit calculations. We thus introduce the perturbed triad^ 

^ p{x)St = {p + 6p{xM (15) 

^In |12l lll^j the triad coefRcient has been caUed p{x) and distinguished from a rescaled p{x). The 
difference is not relevant for most purposes here since our equations are mainly rescaling-invariant. Only 
for a discussion of scale-dependence of correction functions is it necessary to distinguish between these 
variables, as we will discuss later. We thus drop the tilde in order not to overload the notation. 
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which gives rise to a spatial metric of the form qab = \p + Sp\Sab- Here and elsewhere, we 
split off the background part 

p := — p{x)d^x (16) 
Vo J 

with the spatial coordinate volume Vq = J d^x. The latter is assumed to be finite, using a 
compact torus topology of space, but will not appear in the final classical equations. Using 
p, we then define the perturbation 

Sp{x) := p{x) — p such that J 6p{x)d'^x = . (17) 

In the usual notation using the scale factor and the scalar metric mode ip in longitudinal 
gauge we have qab = a^(l — 2ip)6ab, leading to the identification \p\ = and 6p = —2pip. 
Remaining space-time metric components are the lapse function = N + 6N and the shift 
vector A^'*. Their usual notation N = a{l + (f)) gives N = a, 6N = acf) for the lapse function 
in conformal time gauge, while A^*^ is zero for scalar modes in the longitudinal gauge. In 
the absence of anisotropic stress, i.e. non-diagonal components of the spatial part of the 
energy-momentum tensor, (p is not independent of if) but has to agree with it. We will 
make this identification in the final equations, but have to keep N as well as iV^ as free 
Lagrange multipliers in initial steps. 

With a diagonal triad, the connection cannot be diagonal in inhomogeneous situations. 
In fact, for a perturbed triad ()15|) one can compute the spin connection to be 

^.^l^^,^j5p_ ^^^^ 



which is antisymmetric and thus non-diagonal. The diagonal part of A\ is then contributed 
solely by extrinsic curvature, which, again in longitudinal gauge where the shift vector 
vanishes,^ is proportional to a time derivative of the triad and thus diagonal, 

K = k{x)5l={k + 5k{x))5l (19) 

A perturbed connection then has the form 

A\ = ^{k + 5k{x))5l + T\{x) (20) 

split into the perturbed but diagonal extrinsic curvature part and the non-diagonal (in 
fact, antisymmetric) part coming from the perturbed spin connection. The direct calcu- 
lation (fTH|) can easily be seen to solve the Gauss constraint identically. The gravitational 
variables p, k, 6p and 6k are thus constrained only by the diffeomorphism and Hamiltonian 
constraints. 



■^As a gauge choice, this is to be used with care in deriving equations of motion, as described before. In 
this case, one can check that equations of motion of non-diagonal triad terms only express non-diagonal 
extrinsic curvature terms through D^^Niiy Thus, the shift vector can be consistently set to zero. 
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2.3 Constraints 



The gravitational contributions to these constraints in terms of Ashtekar variables are the 
diffeomorphism constraint 

= \ [ d'xN'^Fj^.El (21) 

where F*^ = daA\ — dbA^ + e^^^A->^A^ is the curvature of the Ashtekar connection, and the 
Hamiltonian constraint 

Hg[N] = ^ [ d'xN |det E\-'/' {e.jkF:,E^El - 2(1 + i')KKiEfE'f) . (22) 

Note that, strictly speaking, ()2ip is the so-called vector constraint, i.e. a combination 
of the diffeomorphism constraint which generates spatial diffeomorphisms and the Gauss 
constraint (jl4j) . Since by our choice of variables we are solving the Gauss constraint 
identically, we can use (j2ip as the diffeomorphism constraint. 

Clearly, the matrix symmetry of (fT3j) . (fTHj) and (fT^ will lead to simplifications of the 
Hamiltonian and other constraints. Since constraints are also the places where quantum 
corrections enter, in particular in ()22j) which determines dynamics, one should use simpli- 
fications only after such corrections have been implemented. Specifically, in the effective 
regime, the terms containing inverse powers of the scale factor, e.g. the spin connection 
fll8|) or |deti?|~^/^, will have quantum corrections. If one uses (fT^ first, then the in- 
verse determinant in ()22|) would cancel at the classical level already, thus not acquiring 
any effective corrections upon quantization. Reliable places for correction terms can thus 
only be found in the full expressions. In this paper, we will accordingly insert a function 
of triad components a to take into account possible corrections coming from | deti?|~^/^ 
and a function (3 multiplying spin connection components ()18|) to take into account the 
inverse triad component there. These functions depend on triad components and clas- 
sically a = /5 = 1 in which case we will indeed reproduce the classical equations. We 
will, however, keep a and (3 in our equations without specifying them to demonstrate how 
quantum corrections can easily propagate in more complicated terms when equations of 

*The spin connection is in general more complicated as obtained directly from However, the 

spin connection is not quantized in the full theory because it is not a covariant object and does not 
have physical meaning without gauge fixing. There is thus no tight prescription on how to obtain 
correction terms in perturbative situations. A further difference between our perturbative treatment 
here and the full theory is that we treat extrinsic curvature and the spin connection separately as it 
has been proven useful in homogeneous ^1] and midi-superspace models II 5| . This is not done 
in the full theory where one rather quantizes extrinsic curvature components in the constraint using 

if^ cx I J xe^^^ F^i^-^==^, j yJ\detE\d^ _lg^. Such a treatment is also possible to analyze in 

our situation but would be more complicated. We will later see that corrections from the spin connection 
are much less important than corrections from the inverse determinant explicit in the Hamiltonian con- 
straint. Thus, the precise treatment of extrinsic curvature and the spin connection does not seem crucial 
for phenomenology. 
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Figure 1: Typical behavior of correction functions a, /?, D and a which approach one from 
above for large arguments. For small arguments, the functions are increasing and reach a 
peak value larger than one. 



motion are derived; such functions can be found in The typical behavior is similar to 
the function \p\^^'^d{p) [T7l|TH] used in isotropic models and sketched in Fig. [H For small 
arguments these functions are increasing, starting from a value smaller than one which is 
zero in homogeneous models and for perturbative inhomogenities but can be non-zero if 
non-Abelian features of the full theory and coherent state effects are considered |19| 1201 • 
At an intermediate scale, whose value depends on details of the quantization, they peak at 
a height larger than one and then approach the classical value one from above. The small-p 
behavior is thus strongly modified by non-perturbative effects while the large-p behavior 
is perturbative of the form 

a{p) = 1 + c f with c> and n > (23) 

where ^■p = \/fiG is the Planck length. The values of c and n also depend on the quan- 
tization but are generally positive. This allows the discussion of characteristic qualitative 
effects. It is safest to use perturbation theory above the peaks of correction functions. 

The constraint (|22j) can be simplified, if expressed in terms of the spin connection and 
extrinsic curvature, using A\ = T\ + The first term is then written as 

Fi, = 2d^j:i^ + 2^dyaKi^+ei,,[Ti + ^Ki)[Tl + ^Kt) 

= 2d[an] + maKl] + l^vk {TiK^ + TtKi) + e.,, (F^F^ + ^^KiK^) . (24) 
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Contracting with the triad, we see that the second term 

ei.kdaKlE^El (X e^.kWA « ^^i^^l = 0, (25) 

whereas the cross-term 

vanishes because K is diagonal and F is anti-symmetric. Also, the last curvature term in 
()24|) (quadratic in K), can be combined with the last term of ()22j) yielding 

Hg[N] = - [ d^xNa |det E\-'^' (e,,^^, {pvi) + P^TiVt - KiK^) E^E'J . (26) 

As already announced, here we have included quantum correction functions a and j3 which 
depend on the densitized triad and reproduce classical behavior for a = /3 = 1. 

In the diffeomorphism constraint we have similar simplifications and some terms do not 
contribute after writing out the curvature components explicitly. In our variables, 

Fl,E\ = 2^{p + 5p)d,{k + 5k) + e,,,{Ti + ^Ki)Tl{p + 5p) 
= 27(p + 6p)da{k + 6k) - {Ti + -fKi)djSp 

using ()18|) in the last line. Then, with Ti^djSp = the constraint 

DIN"] = d?x5N\2pda5k - (3kda5p) 

for N"" = results. Note that we did use the expression for the spin connection in the 
second term which could thus receive a correction function (3. Since it is coming from the 
spin connection it should equal /?, but we keep it separate because the diffeomorphism 
constraint is not quantized in infinitesimal form in loop quantum gravity. One rather 
implements finite diffeomorphisms exactly [21] such that no corrections are expected. As 
we will see later, effects of (3 are not important such that the precise prescription here 
does not matter much. Moreover, there is no inverse determinant in the diffeomorphism 
constraint and no need for a correction function a. 

Together with contributions from matter fields, this defines constraints on the basic 
variables. For a scalar field ip with momentum vr and potential V{(p), we have a contribution 

^^[^1 = j d^xN'^ndaip (27) 
to the diffeomorphism constraint and a contribution by the matter Hamiltonian 

H,[N] = [d'xN ( -D^^= + l^E^M^ + ,/\d^\V{^)\ (28) 
W V2 ^^d^ 2 ^/\d^\ ^' ' ^ ' 
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to the Hamiltonian constraint. Again, only the contribution to the Hamiltonian constraint 
contains inverse powers which occur in two different forms. The kinetic part has a single 
inverse determinant which we correct by a function D (as it has been used in isotropic 
models), while the gradient term has additional triad components in the numerator which 
can lead to a different correction function a. Note also that, while formally the combination 
of triad components in the gradient term is the same as that in the gravitational part of the 
Hamiltonian constraint, it is only the symmetric part in a and h which enters the matter 
gradient term but the antisymmetric part in the gravitational constraint. We therefore 
keep the correction functions a and a independent. 

The equations to consider are thus the two constraint equations 

DIN''] = DgIN"] + D^IN"] = (29) 
H[N] = Hg[N] + H^[N] = (30) 

together with the Hamiltonian equations of motion 

f = {f,H[N]+D[N'']} (31) 

for any of the variables p, k, (f, p^ and the fields 6p, Sk, S(f and Sir. The components 
N and N"', i.e. N, 6N and SN°- play the role of Lagrange multipliers for the constraints. 
When computing the Poisson brackets in equations of motion ()31|) one has to keep these 
multipliers as independent variables at this stage, as discussed before. Only N has to be 
specified to fix the time gauge, with the two most common choices N = 1 for proper time 
and N = a for conformal time which we will use here. The fields 6N and 6N"', on the 
other hand, must not yet be fixed to ip or zero, respectively, but be kept independent of 
the canonical fields. 



3 Linear perturbation 

In order to derive linearized equations of motion, we expand the Hamiltonian to second 
order in the field perturbations so as to get linear equations after taking Poisson brackets. 
For the constraint equations themselves, the linear coefficients of 6N and 6N"' will result 
as perturbation equations, accompanied by equations of motion for 6p, 6k, 6ip and 6n, 
which can be combined to two second order differential equations of motion for 6p and 6(p. 
Spatial integrations of terms linear in perturbations give zero because we have split off the 
homogeneous background contributions explicitly in definitions such as (fTTjl for 6p. The 
Poisson structure for these variables can then be computed by inserting them in the full 
action term 

1 f ^3 dAl^^ 3 f 3 dk{x) , , W^dk 3^3 dbk^ 
— \ d^^^" = - / d'-x—^pix) = -^—p + - / d^x—6p 
J dt K J dt n dt K J dt 

such that 

{^,p} = ^ ' {5k{x).5p{y)} = l5{^^y)- (32) 
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3.1 Variations 

All necessary equations follow from variations of the total constraint 

H = Hg[N] + H^\N] + DgW] + D^[iVT 

^YkJ (^-6(A; + 5k) \p + 5p\i + \p + 5p\m 

+ 2 ^^y^'ff, ^ + ^ y d^x^iV^ (2pda5k - mkda5p) + if^[iV] + D^[iV'^] (33) 

by the independent variables iV, 5A^, 5A^", p, 5p, fc, 5ip, vf and ^vr. The background 

shift vector N"- does not appear because it vanishes for homogeneous models, as does the 
background diffeomorphism constraint it would be the Lagrange multiplier of. It would 
have to be considered when perturbing around an inhomogeneous background. 
We then have the following equations: 

= ^ = -^av^P + ™ (34) 

ON K ^ dN ^ ' 

gives the background Friedmann equation which is corrected by higher order terms such 
as / <5?xa5k'^ in the perturbations if higher orders are retained in the expansion. In 
this variation, N is kept independent because it must be varied, but from now on we set 
N = \/\p\ to choose conformal time, derivatives of which will be denoted by a dot (while 
a prime is used for p-derivatives). Then, 

O^'JMl^M UaMk - 3ae (l + 2p^) f + 'l^Hr^ + ™ (35) 

5{5N) K \ V a J 2p p J 5{5N) ^ ' 

gives the first perturbation equation equivalent to (jH}, 

gives the third perturbation equation equivalent to © and 

i . . , (37) 

3Vo op 3Vo dp 

gives the background Raychaudhuri equation. In 

5k = -4^ = ^ (a + a'p) (SNF + 2Nk5k) + t^N¥ (a - a'p - a"p^) (38) 

3 5{5p) ^\p\ \p\-^'^ 

+0 (^5{N 

^ 'a(3i(3 - 2) - APia'p)) VHp + -^WHN + ^-^^-[^^ 



6p3/2 ' > ^" ^ ' 3^ 3 5{5p) 
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which gives the second perturbation equation equivalent to (jH), the term 0{6{N^/\p\)) 
indicates that there are additional terms proportional to 6{N^/\p\) which are not evaluated 
explicitly here. They will cancel exactly in the final equations for the modes used here, 
but would give non-zero contributions if the lapse perturbation and the scalar mode are 
not identified or if other gauges are used. 
Finally, 

relates the connection component k to the time derivative of p or the scale factor a. To- 
gether with the perturbation equation 

5p = —^i—T = 2ap5k + 2akp5p (40) 
6 o{ok) 

which relates the connection component Sk to the time derivative of 6p, it can be used to 
eliminate the extrinsic curvature components. 

For the matter variables we obtain four additional equations, 

7r = --^ = -p'V„{^) (41) 
Vo dip 

which gives the background Klein-Gordon equation, 

5vr = = -p (K,^ {^)5p + V,^^ {^)5^ - a{p)VH^) (42) 

which gives the matter perturbation equation, 

. 1 dH D{p)N _ 

^=v,-d^ = ^rr- (43) 

which relates vf to the time derivative of <^ and 

which relates ^vr to the time derivative of 0. 

Eqs. dnni), pnj)- and (jUj) will be used to eliminate momenta from the equations, 
rewriting some of them as second order differential equations. 



3.2 Metric equations 

We first turn to the more complicated equations obtained by varying with respect to metric 
modes. Here, both the gravitational and the matter part of the constraints contribute, 
whose variations are discussed separately. From now on, we evaluate the variation equations 
only for the case bN"" = (for longitudinal gauge without vector modes) and 6N = 
—6p/2^/\p\ (identifying = ijj). The latter identification implies 6{N^/\p\) = at the 
linearized level which we will use from now on. 
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3.2.1 Gravitational part 

Eq. (j39|l can be rewritten as 

n:=^ = ak (45) 
2p 

where Ti. is tlie conformal Hubble rate. Inserting it in Eqs. (|34p and (|37|) gives the first 
order Friedmann equation 

= (46) 

and the second order Raychaudhuri equation 



\ a J 3Vq dp 



for the background metric. Here, as well as in any equation of motion, the lapse function is 
fixed prior to taking the p-derivative. In other words, any appearance of a time derivative 
implies that a time gauge has been chosen, i.e. the lapse function has been fixed. With 
this in mind, the background Raychaudhuri equation can be written 

Solving Eq. PU]) for 6k we obtain 

6k=-^--k5p, (48) 
zap a 

and inserting k and Sk in terms ofp and 6p in Eqs. (jHUj) and (jHHj) gives the perturbation 
equations 

-—vdp-n ^nil-ap/a)— = — =^-^—— (49) 



a~'d,{-6p + n6p{P + 2a'p/a)) = /^^^ 



(50) 



11 H 

-6p + - {af3{(3 - 1) - 4(3{a'p)) V^6p (1 + 2ap/a) 6p 

a 3 a 

2 



"a^ \a) (2" P + " -4(ap) /a) 5p = Y^~6{6p}~ ^ ' 



for the metric mode Sp. 
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3.2.2 Matter part and energy-momentum 



Rather than using the energy momentum tensor as source, the primary object in a canonical 
analysis is the Hamiltonian combined with the diffeomorphism constraint. The matter 
Hamiltonian is directly related to energy density ^ 

^^"^d^ m ^^^^ 
while contributions to the diffeomorphism constraint give the energy ffux density 

, ' y (53) 

This corresponds to time-time and time-space components of the energy momentum tensor. 
The remaining components, in the absence of anisotropic stress, are pressure components 
which we use here only in the isotropic case. From the thermodynamical definition of 
pressure as P = —bE^jW with energy E^p and volume V , pressure components can then 
be derived through ^ 

A^5v^|detE| 

from the Hamiltonian. 

For the perturbative treatment, we again split these expressions into background and 
perturbation parts such as p and bp. By the chain rule, we have 



P, = -±,^li^ (54) 



where we used 



bN{x) dN J " 6N{x) 6{6N{y))J 

'"^^ ^' Uo dN ^6{6N{x)) Vj "^ysmiy))) ^^^^ 
p^ + 5p^{x) 

for N := Vq~'^ J d^xN(x) and SN(x) = N{x) — N. The last term in vanishes be- 
cause 6H/6{6N{y)) is linear (or in general odd) in perturbations and thus vanishes when 
integrated over space. The remaining terms then define the background energy density 

1 dH^[N] _ ^0 
" W^'Vo dN ~ ' ^^^^ 



'"^Although there is an inverse triad in this equation, it is not quantized as there is no energy density 
operator in an inhomogeneous setting of loop quantum gravity. In any case, such a re-definition would just 
change the relation between energy density and energy-momentum components but not affect the primary 
equations (|49I51|) and lf72|) . 

^In this definition, in contrast with (|47() . the lapse function is independent of the triad. 
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and the linear perturbation 

6p (x) - , 1^1-3/2 ^HAN] _ 

^PA^)- 2\p\^/W, on 5{6N{x))~ ^^^^ 

Thus, 

- = IP^ST^ + IVWlf^Sp . (58) 

Similarly, we obtain 

- _ 1 dP^jN^] 

^'"~7V|p|3/2 dN- " " ^ ^ 

which vanishes for a homogeneous background, and with this 

for the flux. Finally, we have 

^ 3N^\Vo dp " 

and 

for pressure. Note again that the lapse is treated as an independent function at the stage 
of differentiation. Then using 6{N^/\p\) = 0, this gives 

For a scalar field with correction terms in the Hamiltonian, these formulae yield the 
energy-momentum components 

I? = (65) 



2pD'^ \ 3D 
for the background and 

5Tt = (68) 

(69) 



2p2D I 3L) 3 \ D j 3D pD \ 3D 
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for perturbations. 



3.3 Matter equations 

Solving Eq. ()43|) for vf in terms of (p and inserting it into Eq. ()41|) yields the Klein-Gordon 



equation 



V \ D 
for the background scalar field i^. 

Taking a time derivative of Eq. one gets 



^ + P^(l-^^]+pDV,^{ip) = (70) 



p J D p 



where the previous equations (jl^ and (jlHj) have been used. Finally, substituting vf and 
from (j41|) and (j42p . we arrive at the Klein-Gordon equation for the perturbed part of the 
scalar field 

6^ + (^1 - - DaV^6^ + Z^pV^,^^ + (D - D'p)V,^ {^)6p 

3.4 Translation to metric variables 

We can now finally write our equations of motion in familiar form by replacing derivatives 
of the matter Hamiltonian by energy momentum tensor components and by introducing 
the scalar mode = —6p/2p. We keep the variable p rather than expressing it as the scale 
factor squared since this is the basic quantity appearing in our corrections functions from 
quantum gravity. 

The background equations pHjl and (jUj) become 



n"^ = '^app^ (73) 



and 



n = -n' (i - ^) + ^ap (p^ - 3P^) 

using and (jHH). Eq. (jlHI), resulting from (jHHjl . together with (jK^ yields 



(74) 



af3V'ij - -n^ - -n'^ ( 1 - ^ ) = -^p5T° , (75) 



a a \ a J 2 
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Eq. (fFTj). resulting from (jHHj) . together with (jUHj) yields 



\ a J \ 6 a J 6 

+ fl _ 5^^^ + 4 y - 2^) = ^piT: (76) 

\ a \ Ci J a I 2 

and from ()50p together with ()60p we obtain 

da + ^^(2 - /3 - 2aW«)) = -^P-^^a • (77) 
Here, the energy-momentum tensor components are 

= -^^Sip,a (79) 

2pD y 3 D 3\ D J 3 D J pD \ 3 D J ^ 
and is subject to the Klein-Gordon equation 

ip + 2nip - +pDV,^ (v9) = (80) 

for the background and 

+ 2m^ [l - - DoVH^ + DpV,^,, {^)Sip + 2{D - D'p)pV,^ (43)1/' 



for the perturbation, obtained by expressing (j72j) in terms of if). 

Comparison with Eqs. Q, 0, (0), O and (fTUI) shows that the classical equations 
are indeed reproduced when all correction functions equal one. As the classical perturbation 
equations, the corrected ones are scale invariant when the background is fiat. This is 
manifest in the written form since ip is scale invariant and, although the background scale 
factor a? = \p\ appears explicitly, any combination such as a'p where the prime denotes a 
derivative by p is scale invariant, too. Although the correction functions depend on p, the 
derivation shows that they do so only in combinations which are scale invariant, taking 
into account normalizations provided by a quantum state ^3 ^] • 
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4 Covariance 



We have derived corrected perturbation equations in a fixed gauge, which simphfied quan- 
tum and classical calculations. Their space-time covariance is thus not obvious, just as the 
classical equations in the form Q, © and (0) are not manifestly covariant. When classical 
equations are modified by quantum corrections, in particular in a canonical scheme, it is 
not clear whether covariance will be broken. Canonically, space-time covariance is realized 
if the lapse function N and shift vector A^*^ are not restricted by equations of motion but 
can be specified freely as a gauge choice. This is always the case if the Hamiltonian and dif- 
feomorphism constraints form a first class set, i.e. their Poisson brackets among each other 
vanish when the constraints are satisfied. With arbitrary modifications in their terms, this 
is unlikely to remain true, suggesting a breakdown of general covariance. 

The situation in quantum gravity is, however, more general because new quantum de- 
grees of freedom arise, which can absorb some of the restrictions which would otherwise be 
imposed on N and N"". In fact, for an effective description of a canonical quantum theory 
[22] one derives effective constraints such as H^s = (H) and Des = {D) as expectation 
values in suitable states. Since states are described by many, in fact infinitely many, more 
variables (or fields in a field theory) than just the classical ones the effective constraints are 
imposed on all these parameters. Additional variables include, e.g., the spread or deforma- 
tions of wave packets in addition to expectation values identified with classical variables. If 
the quantum constraints preserve the first-class nature of the classical constraints, one has 
by definition {ifefr? D^s} = ^{[H , D]) as a first-class set of constraints. General covariance 
is thus preserved. 

However, effective descriptions not only entail taking expectation values but also a 
truncation of the infinitely many quantum variables to a finite set (which means finitely 
many fields in a field theory). This corresponds, in some sense, to the derivative expansion 
done in effective actions to arrive at a finite sum of local correction terms. In particular, 
in this paper we completely ignored, as a first approximation, all quantum variables and 
correction terms they imply. Such truncations usually lead to effective constraints which 
do not exactly preserve covariance. Nevertheless, in perturbative regimes of quantum cor- 
rections the equations are consistent: One can choose a classically motivated gauge, fixing 
N and A^", and compute corresponding perturbative corrections as we did for the longitu- 
dinal gauge. However, the gauge should not be too special as assumptions could implicitly 
be used which will no longer hold with quantum corrections. This would be the case if 
one used conditions on spatial metric components to achieve a certain gauge, obtained by 
solving gauge transformation equations for lapse and shift. Since gauge transformations 
change themselves when constraints are corrected, such a gauge would not be safe for the 
derivation of corrected equations. Specifying lapse and shift directly, such as A^'^ = in 
longitudinal gauge, is safer because it can be done in the same way with any corrected con- 
straints. This gauge is not complete, but after having computed the corrected constraints 
one can combine the remaining variables to gauge invariant quantities. 

Non-perturbative regimes, such as those below the peaks of our correction functions, 
have to be treated with more care since the classical background geometry is strongly 
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modified there. Usually, additional quantum variables are required to describe the situation 
and to discuss gauge issues fully. (The homogeneous background evolution, on the other 
hand, is safe even in this regime since it is subject to only one constraint. This will 
automatically commute with itself, thus being first class.) 

5 Discussion 

We have presented in detail the Hamiltonian derivation of cosmological perturbation equa- 
tions for scalar modes in longitudinal gauge around a flat isotropic background. The same 
scheme, of course, applies to other gauges and also under inclusion of vector and tensor 
modes and for perturbations around different backgrounds. As in the case of scalar modes, 
due care has to be taken in deciding when a gauge or mode selection is to be specified. 

Since our main interest is to compute corrections from a canonical quantization of 
gravity, typical such correction functions have been included. We have seen how simple 
modifications of the constraints can propagate to more involved corrections of equations 
of motion. We emphasize that we have not presented a complete set of effective equations 
including all possible correction terms. Alternative corrections can arise, and moreover 
gauge issues have to be studied. 

For applications, it is important to note that not only coefficients in evolution equa- 
tions are corrected, but also constraints are modified. Since constraints generate gauge 
transformations, the form of gauge invariant variables changes, too. For instance, it is 
not sufficient to take the classical expression of the gauge invariant curvature perturba- 
tion TZ = ip + T-C6(f/(f and use corrected equations of motion for all variables involved. A 
complete treatment requires correction terms in matter and metric equations as derived 
above, as well as in expressions for the relevant quantities to be related to observations. 
Ignoring any of the ingredients in general can lead to misleading conclusions. Nevertheless, 
some qualitative conclusions can be drawn. For instance, a modified evolution equation 
for TZ can be derived which implies correction terms leading to a slight non- conservation 
of this curvature perturbation [23J. While the quantity TZ itself will have to be corrected 
as the relevant gauge-invariant quantity, implying additional corrections to the evolution 
of curvature perturbations, this is unlikely to happen in such a way that all corrections 
from equations of motion and gauge invariance properties conspire to cancel each other. 
For the precise form of non-conservation, however, all these effects have to be taken into 
account. Still, interesting qualitative effects for cosmological phenomenology have already 
materialized. 

We have started here a program to derive effects systematically and presented a first 
set of corrected constraints as well as evolution equations. A systematic study of different 
gauges and of observable implications is still to be done. The derivation in a Hamiltonian 
formulation as well as the use of Ashtekar variables are crucial for the inclusion of effective 
quantum gravity effects in modified perturbation equations if canonical quantum gravity, 
in particular loop quantum gravity, is employed. 

Primary dynamical objects are then the constraints, rather than Lagrangians, which 
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are modified by quantum effects. Without regarding gravitational parts and all matter 
energy-momentum terms, changes to the classical behavior in an inflationary context have 
been considered in j2ll 123 EE] in a strongly modified regime of background correction 
functions and in j2Z| in a perturbative regime. Ignoring corrections in gravitational parts 
of the equations corresponds to choosing a fiat gauge in which no metric perturbations 
are present. The availability of this gauge choice is based on classical reasoning, and 
has to be reconsidered with gauge transformations generated by the quantum corrected 
constraints. Our treatment is more general since we allowed metric perturbations ip and 
matter perturbations 5ip to be independent before they are to be combined to a quantity 
gauge invariant under the quantum corrected transformations. This allowed us to discuss 
non-conservation of curvature perturbations as a new effect. 

Following the lines of derivations in this paper, basic effects in the constraints then 
translate unambiguously into effects in perturbation equations. Since several of the vari- 
ational equations have to be combined in different manners, even simple modifications 
in the constraints can have complex implications at the level of perturbation equations. 
Modifications one expects on general grounds are regular versions of any inverse power 
of metric variables such as those of p in the Hamiltonian constraint, the spin connection 
and the matter Hamiltonian j2Hl ^1 1211 ^1 higher order corrections as powers of k 
|3n | I3 H l32] and higher derivative terms in space as well as in time j22]. All this gives rise 
to characteristic correction functions which can be computed at least qualitatively. In the 
perturbation equations, coefficients as well as the derivative order of the equations can 
then change and differ considerably from the classical ones in strong quantum regimes. 

We have illustrated this throughout the paper with corrections which are expected 
from inverse power modifications. Those corrections are easiest to implement and to deal 
with because they change only coefficients but not the type of perturbation equations. 
They are also expected to be stronger in inhomogeneous situations [12]. One expects four 
different correction functions, two for the gravitational Hamiltonian and two for the matter 
Hamiltonian. When they equal one, classical behavior is reproduced, while on small scales 
they can differ considerably from one and lead to modified and new coefficients. On very 
small scales, i.e. in regimes where correction functions are not Taylor expandable around 
the value one, cosmological perturbation theory is more difficult to apply. 

There are many effects from quantum gravity in combination, and even different im- 
plementations depending on the quantization scheme used for constraints. An effective 
analysis shows which of the terms are most crucial for physical consequences and should 
be fixed. Other corrections on which the behavior does not depend so sensitively can then 
first be ignored. We can clearly see this from our example, where the correction functions 
j3 from the spin connection and a from the matter gradient term do not play as important 
roles as the functions a and D. This is fortunate in particular for (3 because there is no 
tight prescription for its behavior in the full theory. Also the function /3, which could equal 
(3 or simply one depending on how one deals with the diffeomorphism constraint, only ap- 
pears once in the final equations (j77|) and in a way which does not significantly change the 
behavior given by a-corrections. (For large p, 1 — /5 is negative while —a'p/a is positive. 
Due to the perturbative form ()23|) of the functions as a power series in however, the 
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correction from a is dominant in this regime and determines the sign of the correction. 
For small p, on the other hand, effects from /3 can be more pronounced but perturbation 
theory is more complicated.) The most sensitive corrections at the level of linearized per- 
turbations around flat space are thus those coming from a and D. A phenomenological 
analysis then shows which behavior of these functions is preferred. 
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